The charge shuttle as a nanomechanical ratchet 
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We consider the charge shuttle proposed by Gorehk et al. driven by a time-dependent vohage 
bias. In the case of asymmetric setup, the system behaves as a rachet. For pure AC drive, the 
rectified current shows a complex frequency dependent response characterized by frequency locking 
at fracional values of the external frequency. Due to the non-linear dynamics of the shuttle, the 
rachet effect is present also for very low frequencies. 
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The great burst in the study of Nano-Electro- 
Mechanical (NEMS) devices is unvealing several new per- 
spectives in the realization of nanostructures where the 
charge transport is assisted by the mechanical degrees 
of freedom of the device itself. This is the case, for ex- 
ample, when nanomechanics has been combined with 
single electron tunneling. Important experiments in this 
area are the use of a Single Electron Transistor (SET) 
as dispacement sensor J3] or quantum transport through 
suspended nanotubes [J oscillating molecules 0, 0, S 3 
and islands lal . On the theoretical side, several works 9|, 
ESIIlIillllllllllllllllliillllii have already 
highlighted various aspects of the role of mechanical mo- 
tion on single electron tunneling. 

An exciting prototype example of mechanical assisted 
SET device has been proposed by Gorelik et al. and 
named Single Electron Shuttle. The authors of Ref. 
predicted that a SET with an oscillating central island 
can shuttle electrons between the electrodes leading to 
low noise transport [lol . lrRll^ . Although the realization 
of the charge shuttle is difficult experimentally, promising 
systems are molecules in break junctions 0, 0, or 
silicon structures 0, 0| . 

One of the advantage of this self-oscillating structure 
is the fact that one can generate very high frequency me- 
chanical oscillation with static voltages. The fact that 
the system has an intrinsic and stable oscillating mode 
as the result of a static voltage suggests that the applica- 
tion of an oscillating voltage may lead to new interesting 
effects, related to the interplay between the external AC 
drive and the internal frequency of the device. Moreover, 
as the non-linearities of the dynamics have an important 
role, this interplay should emerge in a wide interval of 
the ratio of the two frequencies. Aim of this Letter is 
to study a shuttle driven by a time-dependent applied 
bias. The most interesting situation is when the system 
is asymmetric. We will show that the structure acts as 
a ratchet in which the forcing potential is generated 
in a self-consistent way. We study in some details the 
case in which the external bias is AC and report a quite 
rich behaviour as a function of applied frequency. We 
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FIG. 1: Left panel: Schematic of a charge shuttle. Right 
panel: energy diagram for the SET. The diagonal lines indi- 
cate the thresholds for the vanishing of the four rates Vfl, 
Ffk, Ftl, and Ttr- The two dots indicate the state of the 
system during shuttling at fixed voltage V. 



find clear indications of frequency locking. The result- 
ing DC current may have both signs, depending on the 
value of the frequency. The response to an AC field can 
give strong indications on the motion of the central is- 
land, even when the system is very far from the shuttling 
instability. A sizeable ratchet effect is present down to 
frequencies much smaller than the mechanical resonat- 
ing frequency, due to the adiabatic change of the equilib- 
rium position of the grain. In a very recent experiment, 
Scheible and Blick |2l| . already observed similar results 
to those presented in our work. 

The single electron shuttle, shown in Fig. ^ is a SET 
where the central island can oscillate between the two 
leads 01 ■ The central island is subjected to an elastic re- 
coil force, a damping force due to the dissipative medium, 
and an electric force due to the applied bias. The is- 
land is connected to the left and right leads through tun- 
nel junctions with resistances Rl{x) = RL{0)e^^^ and 
Rr{x) = RR{0)e~^^^, (A is the tunnelling length and x 
is the displacement from the equilibrium position in ab- 
sence of any external drive). The capacitances Cg and 
C couple the island to the gate and to the two leads, 
respectively. We neglect small variations of the capaci- 
tances due to the motion of the island as this dependence 
is weak when compared with the exponential dependence 
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of the tunneling resistances on x. The system is biased 
symmetrically at a voltage V{t) (Vr = —Vl = V/2), the 
charging energy is Ec = e^/2Cs where e is the electron 
charge, = Cg + 2C, and the gate charge Qg is CgVg. 
For simplicity, we consider the case of low temperatures 
(T <C Ec), charge degeneracy {Qg = e/2), and voltages 
\V\ < Ec/e. In this regime the grain can accommodate 
only n = or 1 additional electrons (see right panel of 
Fig.^. All the properties of the nanomechanical ratchet 
are captured already at this level. 

In the simplest approximation the dynamics of the cen- 
tral island is described by Newton's law 



x{t) = -ujI x{t) - J x{t) + 
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n{t) 



(1) 



of the grain, loo is the oscil- 
7 is a damping coefficient and 



Here m is the mass 
lator eigenfrequency, 
L is the distance between the two leads. In the 
regime of incoherent transport the (stochastic) evolu- 
tion of the charge —en{t) is governed by the follow- 
ing four rates [23]: TpL = \eV{t)/4Ec\TL{x)QiV) and 
= \eV{t)/4:Ec\TR{x)e{~V) for n = ^ 1 tran- 
sitions, Ftl = \eV{t)/4:Ec\TLix)ei-V) and Ttr = 
\eV{t)/4:Ec\TR{x)e{V) for n = 1 ^ transitions. 
Here FL, FR, TL, and TR stands for From/To and 
Left/Right, indicating the direction for the electron 
tunneling associated to the corresponding, Fi//j(a;) = 
[R[^/r{x)C]~^ and Q{t) is the Heaviside function. The 
current / is then determinated by counting the net num- 
ber of electrons that have passed through the system in 
a given time t. 

In most of the paper we consider the case of an os- 
cillating voltage bias V{t) — VoSin(wt) where the 
interplay of the two frequencies lo and cuo is crucial. If 
the system is perfectly symmetric no direct current can 
be generated, since this would break the left /right sym- 
metry. We then concentrate on the asymmetric case 
Rft.{0) / Rl{0) ^ 1. We performed simulation of the 
stochastic process governed by the four rates defined 
above and by Eq. The results presented are ob- 

tained by simulating 10^ events of tunneling for each 
plotted point. After a transient time the system reaches 
a stationary behavior. In Fig. |21 we show the station- 
ary DC current as a function of the frequency of the 
external bias. The rich structure shown in the figure is 
generic, we observed a qualitatively identical behaviour 
in a wide range of parameters. The existence of a di- 
rect current as a result of an applied periodic modulation 
shows that the charge shuttle, whose stochastic dynamics 
has been defined above, behaves as a ratchet |22|. Since 
our system is non linear, the external driving will affect 
the dynamics also for values of co very different from the 
natural frequency tUo- Note that in this model the non- 
linearities are intrinsic to the shuttle mechanism. They 
are not due to a non-linear mechanical force, but they 
stems from the time dependence of n{t). As it is evident 




FIG. 2: Current as a function of the frequency for e = 0.5, 
j/lOo = 0.05, F/lJo = 1, and Rr/Rl = 10. The result of the 
simulation of the stochastic dynamics (points) is compared 
with the approximate la (full line). In the small frequency 
region, enlarged in the inset, several resonances at fractional 
values of u>o appear. We also show (dotted) the analytical 
result from Eq. Q in the adiabatic limit. The triangular dot 
indicates the numerical solution of the adiabatic equations. 
Lower inset: current noise from the simulation (points) and 
analytic result (dashed) for the static SET. 



from Fig. |21 the ratchet behaviour is present also in the 
adiabatic limit lu/lOq I. In addition a series of reso- 
nances, due to frequency locking |25| when lo ^ Uo q/p, 
with q and p integers. In this case the motion of the 
shuttle and the oscillating source become synchronized 
in such a way that every q periods of the oscillating 
field the shuttle performs p oscillations. In Fig. |21 we 
also report the low frequency current noise (lower inset). 
The net DC current results from large cancellations be- 
tween positive and negative contributions. But the cur- 
rent noise is always positive, thus it docs not cancel. We 
find that the noise remains very close to the value for a 
SET S = e^TLTR{Tl + Tr)^/{Tl + Tr)^ shown with a 
dashed line in Fig. |21 At resonance, more ordered trans- 
port is realized and current becomes less noisy [T^ . In 
the following we provide a more detailed analysis by solv- 
ing the problem analytically in some tractable limits and 
by analyzing in more details the behaviour of the shuttle 
in the frequency locked case. 

If the electric force is much smaller than the mechanical 
one, e = eVo/ {ujginXL) <C 1, one can take into account 
the force generated by the stochastic variable n{t) only on 
average by substituting into Eq. ^ its mean: (n(t)) = 
P{t) where P{t) is the probability to have occupation 
n = 1 in the grain. The charge dynamics in the central 
island is then described by a simple master equation |2^ 



p(i) = -ri(t)p(i) + F2(t), 



(2) 



with Fi(t) - \eV{t)/AEc\{TL[x{t]+TR[xit)]) and 
F2(i) = \eVit)/AEc\iTL[x{t]Q[vit)]+TR[xit)]Q[vit)]). 
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The rates depend on the time t through the voltage V{t) 
and the position of the grain x{t). The instantaneous 
(average) current through the structure is 

Iait)/e={l-Pit))TFL-P{t)rTLit) , (3) 

where the subscript in the current indicates that the fluc- 
tuations of the force acting on the shuttle, due to the dis- 
crete nature of the charge tunneling, arc neglected. As 
shown in Ref.^the shuttle instability, at constant bias, is 
controlled by the ratio eWo/7, we can thus assume both 
parameters small, but their ratio arbitrary. 

In the adiabatic limit (lu <ti lUo and ecuo/j <C 1) it 
is possible to find an approximate solution of Eqs. 
and In this case the position of the grain is 

given by the local stationary solution of Eq. x/A = 
eT2{x, V)/Ti{x, V). Solving this equation in lowest order 
in e one obtains 

Lyuj ^ ijJo) = e o ; 7^ — ■ 4) 

The corresponding value is shown in Fig. [3 The (small) 
difference with the full numerics is due to the expansion 
in e. By solving numerically the equation for the local 
equilibrium position of the grain, one obtains the result 
indicated by the triangular dot at = in the inset 
of Fig. El [2y|. If £7/(^0 is larger than the critical value 
for shuttling, the behavior is completely different. The 
grain will oscillate at its natural frequency lOq with the 
amplitude slowly modulated by V{t). The modulation 
will be small, since for small e the effect of a change in 
the value of V affects the mechanical motion only after 
several oscillations. In this case the rectified current is 
much smaller than in the adiabatic limit. 

We now discuss in more details the dependence of the 
current on the external frequency. The most prominent 
structure, observed also in the experiments of Rcf. i2ll is 
present in our simulations at w ~ cJo, and it corresponds 
to the main mechanical resonance. The current changes 
sign accross the resonance. This behaviour is due to the 
phase relation between the driving voltage and the dis- 
pacement of the grain. We verified this conjecture by 
solving Eq. ||2Jl with x{t) = Asin(a;t — (p). By calculating 
the current as a function of (p and using the usual res- 
onant dependence for 0(w) = arctan[(tj^ — uj'^)/ujrj\ we 
could reproduce qualitatively the behavior of Fig. [21 

From Fig.[21it is clear that additional structures appear 
also for Lj « lOoIq (magnified in the upper inset), with 
(7 = 2,3,... (the numerical results indicate that, except 
for the fundamental frequency, even q are favorite with 
respect to odd ones). As we already anticipated, the 
motion of the shuttle and the oscillating source become 
synchronized at commensurate frequencies whose ratio is 
pI q. This ratio, known also as the winding number, can 
be defined more precisely as 

w = lim 6{t)lujt (5) 




FIG. 3: The top panel shows the winding number obtained 
from the stochastic simulation (light points) and from the 
average approximation (black points). In the bottom panel 
Acj) is obtained from the average approximation (dashed line), 
from the stochastic simulation using w from Eq. (light 
points). The variance has also been calculated from the sim- 
ulation in the case of u; = 1 and 10 = 2 (long dashed line). 
The parameters are the same as in Fig.|21 

where 9{t) is the accumulated angle of rotation of the 
representative point (i,i) [0{t)/2'K gives the number of 
oscillations performed by the shuttle during the time 
t]. When the system is frequency locked at a wind- 
ing number w, it is possible to define the phase shift 
4>{t) = Oit) — wujt 123. After a transient time for perfect 
locking (j) should not depend on t (apart from a small 
fluctuation if the motion is not perfectly harmonical). 
An additional important quantity to analyze is thus the 
phase shift variance (Ac/))^ = (i^^) — (c/>)^. This is calcu- 
lated by sampling 20 points per cycle over 10^ cycles. For 
a given w the smaller is the value of A0, the better the 
system locks to that external frequency. The numerical 
results for w, <j), and are shown in Fig. Oand Fig.^ 

In Fig. 121 we show the dependence of the winding num- 
ber as a function of the external frequency (top panel) 
together with the analysis of Ac/) (lower panel) calcu- 
lated from the stochastic simulation and from the av- 
erage approximation (short dashed line). The locking at 
rational winding numbers is confirmed by the presence 
of plateaux of decreasing width. As expected the most 
stable plateau is at w = 1: the system locks very well at 
this frequency. We find that this holds up to very high 
frequency in the average approximation, where w = 1 
seems the only possible winding number. The stochastic 
simulation would indicate instead that for uj/lJo 1-3 
the locking with w = 1, as defined from Eq. (jsj, is no 
more established. However, by studying Ac/) for w = \ 
in the whole frequency range, one actually obtains that 
a correlation is always present [i.e. Ac/) < 7r/-\/3 ~ 1.81) 
even when Eq. jSJ gives a value of w different from one. 
The stochastic fluctuations thus unlock the shuttle glob- 
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FIG. 4: Phase shift iji (same parameters and same notation of 
Fig-HJ- The phase shift is shown only when A(^(ti;) < 1.7 (cfr. 
Fig. 13 . When locking is achieved in presence of the stochastic 
fluctuations <j} remains the same of the average simulation. 

ally, but not locally. Looking at the presence of local 
locking at other winding numbers we find that for in- 
stance w = 2 is clearly present for lu > 1 and reversely 
w = 1 is present around to — 1/2 (sec long dashed lines 
in lower panel of Fig. For global locking, only one 
phase variance is minimal. It corresponds to the "dom- 
inant" winding number. The presence of correlations of 
other winding numbers may indicate partial locking at 
these winding numbers (as in the region 0.6 ^ w ^ 0.9 
for w = 1 and 2) or the contribution of higher armonics 
of x{t) (as in the region uj > 1). 

The dependence of the phase shift, see Fig. 21 around 
each resonance is also remarkable. It is very similar to 
the behavior of a forced harmonic oscillator, but instead 
of evolving from to tt it goes from 0o to (/)o-l-7r, where 4>o 
depends on the resonance. With the aim of understand- 
ing the behaviour of the current close to the resonance, 
we also considered the adiabatic limit, and looked for 
harmonic oscillations superimposed to the adiabatic so- 
lution given before. This can be verified numerically by 
searching a periodic solution of period 27r/a; for P{t). We 
find that in this case the direct current is non vanishing 
and that it depends strongly on (j). Using the phase de- 
pendence given by the full numerical calculation we could 
reproduce the shapes of the resonances. 

In this Letter we showed that the charge shuttle un- 
der time-dependent driving behaves as a ratchet. In the 
case of a AC bias the response of the shuttle is rather 
rich due to the non-linear dynamics of the grain. All 
the results obtained here can be tested experimentally, 
indeed the very recent paper by Scheible and Blick |21| 
already reports on some of the properties related to the 
main resonance. 
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